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INTRODUCTION 


In  conservative  vibration  and  static  buckling  problems,  a  variety 
of  methods  of  approximation  with  their  theoretical  foundations  residing 
in  the  calculus  of  variations  have  been  applied  to  determine  numerical 
values  of  natural  frequencies  of  vibration  and  buckling  loads.  The 
familiar  methods  of  Ritz  and  Galerkin  are  probably  the  most  prominent, 
and  the  theoretical  basis  for  the  application  of  these  methods  to  self- 
adjoint  boundary  value  problems  is  rather  well  established.  The  prin¬ 
ciple  of  stationary  complementary  energy  has  also  been  used  for  such 
problems,  and  despite  the  accuracy  of  the  method,  which  is  particularly 
advantageous  in  computing  the  lowest  as  well  as  higher  eigenvalues,  the 
technique  has  not  received  widespread  attention,  possibly  because  of  the 
extra  effort  required  in  computing  expressions  for  stresses  (bending 
moments  in  beam  and  plate  problems)  once  an  apposite  set  of  coordinate 
functions  has  been  selected. 

For  nonconservative  problems  of  elastic  stability,  in  which 
nonself-adjoint  boundary  value  problems  arise,  the  question  regarding 
convergence  of  the  Ritz  and  Galerkin  procedures  has  been  answered  in 
the  affirmative  only  for  special,  albeit  fairly  broad  classes  of 
problems.  Indeed,  the  method  of  complementary  energy  has  received  no 
attention  at  all.  This  seems  rather  unfortunate  since  in  nonconserva¬ 
tive  stability  problems  an  accurate  knowledge  of  the  variation  of 
frequency  in  at  least  the  first  two  modes  with  the  load  parameter  is 
essential  for  the  effective  detennination  of  the  value  of  the  critical 


load.  In  this  report  we  shall  apply  to  three  nonself-adjoint  boundary 
value  problems  these  various  techniques  in  conjunction  with  a  formal 
variational  expression  and  the  so-called  adjoint  variational  principle, 
and  their  relative  merits  will  be  compared  and  contrasted. 


BECK'S  PROBLEM 


Let  us  consider  perhaps  the  best  known  problem  of  the  theory  of 
stability  of  non-conservatively  loaded  structures,  namely.  Beck's 
problem  [1].  This  problem  consists  of  determininn  tho  smallest  load  P 
that  is  applied  to  the  free  end  of  a  cantilever  beam  and  which  remains 
tangent  to  the  deformed  axis  of  the  beam  throughout  the  process  of 
deformation,  such  that  the  amplitude  of  the  resulting  oscillations  of 
the  beam  increases  exponentially  with  time.  Loss  of  stability  in  this 
way  is  often  called  flutter.  The  differential  equation  and  boundary 
conditions  of  Beck's  problem  are 


^  9^w  +  oA  w  =  0,  o<x<ii, 

4  ^ 

9r  8x  at 

w(o,t)  =  ^  (o,t)  =  a__w  (?.,t)  =  a  w  (Ji,t)  =  0, 

-  -2' 

ax  ax  ax 

where  wC",”)  is  the  transverse  deflection  of  the  beam  at  the  point 
X  io^^)  and  at  Lime  F  (F>o) ,  El  is  the  flexural  rrgidity,  o  the 
density,  A  the  area  of  the  cross  section,  and  i.  the  length  of  the 
beam. 


CD 


(2) 
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It  is  convenient  to  introduce  dimensionless  time  and  space 


coordinates.  Thus,  we  set 

X  =  ^x,  t  =  ct, 


so  that  (1)  and  (2)  may  be  written  as 


3^  3^  3^ 

3x*  3x^  3t^ 

w(o,t)  =  w' (o,t) 


=  0, 

=  w‘’(l,t) 


0<X<1, 


0. 


(3) 

(4) 


where  the  prime  symbol  '  denotes  differentiation  with  respect  to  x. 
Next  we  asstune  that 

w(x,t)  =  y(x)e^“^. 


(5) 


where  w  tf.e  natural  frequency  parameter.  Substitution  of  (5)  into 
(3}  ^nd  (.4)  yields  the  following  nonself-adjoint  boundaiy  value 
problem: 


Tx: 

y  (X) 


y  Cx)  - 


uV(x)  =  0, 


0<X<1, 


(6) 


y(o)  =  y'(0)  =  y”(l)  =  y”’(l)  =  0. 


(7) 


It  is  our  objective  to  determine  from  (6)  and  (7)  the  variation  of 
the  frequency  parameter  u  as  a  function  of  the  load  parameter  Q.  Since 
the  eigenvalues,  i.e.,  the  natural  frequencies,  of  nonself-adjoint 
boundary  value  problems  may  be  complex  numbers  (recall  that  self-adjoint 
boundary  value  problems  c?4i  have  only  real  eigenvalues),  the  phenomenon 
of  flutter  can  occur.  This  means  that  if  the  imaginary  part  of  w  is 
negative,  then  the  amplitude  of  the  oscillation  will  grow  exponentially 


in  time.  Generally,  in  those  nroblems  in  which  the  kinetic  method 
must  be  applied,  the  frequencies  of  the  various  modes  are  real  for 
sufficiently  small  values  of  Q  and  vary  as  the  magnitude  of  the  load 
increases,  the  frequency  of  the  first  mode  eventually  increasing  and 
that  of  the  second  mode  decreasing  until  the  two  coalesce  at  a  value 
of  Q  denoted  by  For  Q'^Qcr*  “  is  a  real  number,  but,  for 

is  a  complex  number  and  flutter  occurs,  Ber,k  [1]  solved  (6)  and 
(7)  exactly  and  found  that  the  value  of  the  critical  load,  is 

Qcr  =  2.031ir2. 


C8) 


It  shall  be  our  intention  here  to  apply  several  methods  of  approxima¬ 
tion  to  Beck's  problem  in  order  to  gain  some  understanding  regarding 
the  accuracy  of  the  approximations,  having  available  the  exact  value 
in  (8)  to  serve  as  the  basis  for  comparison.  Ultimately j  we  will 
apply  those  methods  which  appear  to  be  more  accurate  or  most  convenient 
to  two  problems  for  which  the  exact  values  of  the  critical  loads  are 
not  known. 


THE  RITZ  METHOD 

The  methods  of  approximation  are  to  be  based  upon  certain  varia¬ 
tional  formulations  of  the  fundamental  problem  in  (6)  and  (7),  Our 
first  approach  will  parallel  variational  formulations  used  by  Leipholz 
[2]  and  Mote  [3],  We  mjiltiply  (6)  by  <Sy(x)  and  form  the  integral  of 
the  result: 

/^[y’^^Cx)  <•  Q  y’'(x)  -  uV(x)]'5y(x)dx  0. 

0 
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Integrating  by  parts  and  using  the  boundary  conditions  (7),  we  obtain 

<5  i -  Q(y’)^  -  “V]dx  +  Qy'(l)6y(l)  0. 

2  o  (9) 

The  nonconservative  nature  of  the  load  is  manifested  in  the  presence  of 

the  variational  term  Qy*  (1)  ‘Sy(i)  in  (9),  A  complete  functional  in  the 

usual  sense  does  not  exist  here  because  the  nonconservative  load  Q  does 

not  possess  a  potential.  Relations  of  the  type  shov/n  in  (9)  are  often 

called  "variational  principles."  '  lis  usage  is  not  consistent  with 

that  of  the  classical  calculus  of  variations  because  there  is  no 

functional  wnlch  is  stationary. 

Leipholz  [4]  -  [8]  has  demonstrated  that  Galerkin’s  procedure  in 
various  forms  can  be  apnlied  to  stability  problems  of  the  type  under 
consideration  here.  We  shall  assume  that 

N 

y(x)  =  y(x)  =  I  y  (x), 

n=l  (10) 

where  the  a^^'s  aro  unknown  constants  and  the  coordinate  functions  y^Cx) 

are  selected  so  as  to  satisfy  all  the  boundary  conditions  in  (7),  In 

particular,  suppose  that  we  choose 
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Substitutjor.  of  (11)  into  (12)  yields 

n(n+l)A^  >  (n-*-lXn+2)Bj^  =  -  (n+2)  (n+3)Cj^ 

(n-l)nAri  ♦  n(n+2)Bjj  =  - (n+2)  (n+3)C„, 

from  which  we  find 


(n+3)(n+2)  ^ 
■  n(n4l) 


=  -2 


(n^3) 
(n+1) 


C„. 


It  is  convenient  to  set  =  n(n+l)/2,  so  that  (11)  becomes 


(13) 


yn(x)  =  i  [  (n+2)  (n'^3)x”'*^^-2n(n'*-3)x’'^^-«n(n+l)x’'*  ],  nM. 

Let  us  now  consider  (9)  and  the  sequence  of  steps  followed 
formally  in  the  application  of  the  Ritz  method.  If  we  substitute 
^10)  idto  (9) ,  we  find 

I  I  %(/  (yaVK*-  Q  Xn  Xk  -  *  Q  ynCDyk^^)^-^^!.  = 

n=l  k=l  0 

but  since  the  6aj.  are  arbitrary  we  arrive  at 


I  ^l\k  *  Q  “  c  r]  =  0, 

n=l 


(14) 


where 


\k=yn(l)yk^^)-/'  yn(^)ykW^x, 

0  o 

Cnk=  /Vn(-’^)yj,W^^- 
0 

Now  (14)  represents  a  system  of  N  algebraic  equations  in  the  N  unknowns, 
aj^,  and  the  condition  for  the  existence  of  a  nontrivial  solution  is  that 


1  •  .  • 


(let  (A^j.  >  Q  =  0. 


‘^nk 


nk^ 


(15) 
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(15)  was  solved  numerically  for  N  =  2.  The  value  of  the 
critical  load  parc'meter  was  found  to  be 

Q^.j.  =  2.19311% 

which  is  to  be  compared  with  the  exact  value  in  (8). 


(16) 


THh  GALORKIN  PROCEDURE 


We  may  carry  the  analysis  forward  a  little  more  to  show  that  for 
(9)  and  (13)  the  Ritz  and  Galerkin  methods  are  identical.  In  view  of 
the  fact  that  the  yj^(x)'s  given  in  (13)  satisfy  the  boundary  conditions 
(7),  we  may  integrate  by  parts  to  obtain 

^lk  =  /Vn' (x)yj' (x)dx  =  /Vn 
0  0 

^nk  “  “  / Vn(x)yk(x)dx  =  /V^' (x)yi^(x)dx. 

0  0 

Substitution  of  these  last  two  results  into  (14)  leads  to 

I  %  /^y^  (X)  +  Q  yn'(x)  -  w^yn(x)]y,^(x)dx  =  0, 

n=l  0 

or,  by  virtue  of  (10), 

[5'®(x)  ♦Qy"(x)  -  (*>^y(x)]y.  (x)dx  =  0, 

0  (17) 

which  we  recognize  as  the  method  of  Galerkin.  Eq.  (17)  leads  to 

precisely  the  same  numerical  result  as  did  the  Ritz  approach  as 

embodied  in  (14).  In  practice  it  is  generally  easier  from  the 

point  of  view  of  the  numerical  work  involved  to  deal  with  (17) 

rather  than  (9). 


Liivr.'soN’s  .'itniuD 


Levinson  [0]  obtained  an  approximate  solution  for  Bock's  problem 
by  modifying  the  fo.nn  of  (0)  and  then  making  the  numerical  calculation, 
ive  shall  record  here  only  the  outline  of  his  technique. 

The  identity 

y'(136y(l)  =  5[y'(l)y(l)]  -  y(l)6y'(i) 
is  inserted  into  (0),  the  result  being 


/^(Cy")^  -  Q(y')"  -  a,V]dx  +  Qy'{l)y(l)}-  Qy{l)6y'(l)=0. 
^  0 


It  was  subsequently  argued  that  this  last  equation  is  equivalent  to 


/^[(y")^  “  Q(y')^  -  +  Qy‘(i)y(l)  =  o 


with  the  auxiliary  condition 


5y  (1)  =  0. 


(18) 


(19) 


Let  us  assume  that  y(x)  is  again  approximated  by  (10)  and  (13), 
Substitution  of  (10)  into  (18)  leads  to 
N  N 

'2 


n=l  k=l  0 


N  N 
n=l  K=1 


or 


N  N 


J,  “n  t'A'K'  ■  Kh  -  "^Vki'**  * 

n=l  K=l  0 

N  N  N  N 


n=l  k=^l 


n=l 


k=l 


(20) 
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But  the  auxiliary  condition  (19)  leads  to 

N 

I  y^CD  =  0, 

n=l 

so  (20)  reduces  to 

^  Q[y‘(l)yk(l)  -  /Vny>]  - 

n=l  k-1  0  0 

o 

which  is  nothing  more  than  (14), 

Therefore,  Levinson's  method  leads  to  exactly  the  same  frequency 
equation  as  do  the  Ritz  and  Galerkin  procedures.  In  general,  it  would 
seem  that  the  Galerkin  approach  would  be  the  easiest  to  apply,  provided 
that  the  admissible  functions  satisfy  all  the  geometric  and  natural 
boundary  conditions  of  the  given  problem.  Leipholz  [2]  has  discussed 
how  the  Galerkin  procedure  must  be  modified  in  the  event  that  some  of 
the  geometric  or  natural  boundary  conditions  are  violated  by  the  form 
of  the  selected  coordinate  functions  y^C^)* 

THE  METHOD  OF  COMPLENfENTARY  ENERGY 

Washizu  [in]  -  [12]  has  applied  the  principle  of  stationary 
complementary  energy  to  some  conservative  vibration  and  stability 
problems.  This  technique  was  also  discussed  by  Reissner  [13],  Libby 
and  Sauer  [14],  and  van  de  Vooren  and  Greidanus  [15],  When  applied 
to  conservative  problems,  the  method  based  on  the  principle  of 
stationary  complementary  energy,  according  to  [12]  and  [15],  is 
really  identical  to  Grammel's  method  [16],  which,  in  theory,  involves 
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an  integral  equation  formulation  of  the  eigenvalue  problem  and  the 
fireen’s  function  associated  with  the  differential  operator  and 
boundary  conditions  of  the  original  problem.  Leipholz  [17]  -  [18]  has 
discussed  the  convergence  of  Grammel's  procedure  and  has  presented  an 
extended  version  of  it.  Using  a  two-tenr  approximation  in  certain  free 
vibration  problems  associated  with  the  flexural  deformation  of  a  beam, 
Washizu  showed  that  the  principle  of  stationary  complementary  energy 
led  to  an  excellent  approximate  value  for  the  lowest  frequency  and  an 
approximate  value  for  the  frequency  of  the  second  mode  that  represented 
a  substantial  improvement  relative  to  the  corresponding  value  obtained 
by  means  of  the  familiar  Ritz  procedure.  Consequently,  since  accurate 
values  of  both  the  first  and  second  modes  are  required  for  the  deter¬ 
mination  of  the  critical  load  parameter  in  those  nonconservative  sta¬ 
bility  problems  in  which  the  kinetic  method  must  be  used,  it  seems 
reasonable  to  speculate  that  the  method  described  by  Washizu  could  be 
a  useful  tool  for  obtaining  an  approximate  solution  of  such  problems. 

In  order  to  bring  into  play  the  concept  of  complementary  energy, 
we  set 

M(x)  *  -  y"(x), 

where  M(x)  is  related  to  the  bending  moment,  so  that  (9)  becomes 

4  -  Q(y')2  >  a,V]dx  +  Qy’(l)6y(l)  =  0, 

^  o  (21) 

under  the  subsidiary  conditions 

m"(x)  -  Qy^\x)  *  &)^y(x)  =  0,  o<x<l,  (22) 

y(o)  =  y’(o)  =  M(l)  =  m' (1)  =  0.  (23) 
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Note  that  (22)  and  (23)  are  equivalent  to  (6)  and  (7),  In  (21),  we 
shall  consider  the  variation  of  both  M(x)  and  y(x). 

For  y(x)  we  once  again  intend  to  use  (10)  and  (13),  and  the 
expression  for  M(x)  shall  be  obtained  by  integrating  (22)  subject  to 
the  last  two  conditions  in  (23).  Now  substitution  of  (10)  into  (22) 
leads  to 

•  I  N 

M  Cx)  =  I  a^M^'(x), 
n=l 

where  we  define 

M^'(X)  =  Qy;'(x)  -  a,2y^(x),  M„(l)  =  M^^d)  = 

In  view  of  M(l)  =  m'(1)  =  0,  two  integrations  of  (24)  yield 

N 

M(x)  =  I  a^MnCx), 
n=l 

whereas  two  integrations  of  (25)  lead  to 

Mn(x)  =  Ql/nW  ♦  yn(l)(l-x)  -  y^(l)]-w^  Ax-t)yj^(t)dt 


0. 


(24) 


(25) 


(26) 


=  Q  fn(x)  -  ‘»^gn(x). 


(27) 


where 

f^(x)  =  n-(n*3)x 


[(n*2) (n+3)-2n(n^3)x  +  n(n+l)x^]. 


2(3n>10) 

®n^^^  '  (n>4)  (n+5) 


6x  x"*^  / 
n+4  *  2  \  ” 


2nx  n(n+l)x 


n'*-4  (n+4)  (n+5) 


1) 

S)l  • 


Upon  inserting  (10)  and  (26)  into  (21),  we  obtain 
N  N 

^2"  I  I  ®n®k  /^[M  (x)Mj^(x)  -  Qy’(x)y*(x)  -  wVn(x)yj^(x)]dx  + 
n=l  k^l  0 

N  N 

+  Q  I  I 

n=l  k=l  " 
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I’erfcnning  the  indicated  variation  in  this  last  result,  we  find 

N  j 

L  a„{  /  [Mn(x)Mj.(x)  -  Qy;(x)y^(x)  -  a.2yn(x)y^(x)]dx  ^ 
n=l  o 


*  Qy^CDy^d)}  =  0. 


The  numerical  solution  of  the  above  equation  for  N  =  2  leads  to 


=  2.095n 


for  the  critical  load  in  Beck's  problem.  With  respect  to  the  Ritz 
and  Galerkin  methods,  the  complementary  energy  method  trades  increased 
complexity  for  possibly  better  results* 


Till;  METHOD  OF  THE  ADJOINT  VARIATIONAL  PRINCIPLE 


In  the  preceding  paragraphs  we  applied  the  methods  of  Ritz, 
Galerkin,  and  complementary  energy  to  a  nonself-adjoint  boundary  value 
problem  by  means  of  the  variational  expression  given  in  (9).  The 
numerical  results  for  the  value  of  the  critical  load  obtained  in  this 
fashion  with  a  two-term  approximation  were  slightly  greater  than  the 
exact  value  given  in  (8),  the  best,  approximation  being  given  by  the 
method  of  complementary  energy. 

Since  the  applied  load  in  Beck's  problem  is  nonconservative,  it 
does  not  possess  a  potential,  cind  thus  no  functional,  viz.  no  Lagrangian 
function,  exists.  Consequently,  a  complete  variational  formulation  in 
the  classical  sense  cannot  be  stated  for  this  problem.  Incomplete 
variational  formulations  of  the  type  appearing  in  (9),  while  they  often 
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provide  useful  numerical  schemes,  as  we  have  seen  in  our  calculations 
above,  leave  something  to  be  desired  in  the  eyes  of  the  mathematical 
purist  who  prefers  a  sense  of  completeness  in  the  theory  underlying 
his  numerical  operations.  In  this  regard,  the  papers  by  Finlayson  and 
Scriven  [19]  and  [20]  are  well  worth  reading.  Since  no  functional 
exists  for  Beck*s  problem,  m  are  not  finding  a  stationary  value  in 
(9).  Jiq.  (9)  serves  merely  as  a  systematic  means  for  acquiring 
approximate  values  for  the  critical  load.  As  we  demonstrated  above, 
(9)  is  entirely  equivalent  to  Galoi'kin’s  method,  which  is  perhaps  the 
best  known  example  of  a  class  of  techniques  known  as  the  method  of 
weighted  residuals. 

in  the  following  paragraphs  we  shall  introduce  a  boundary  value 
problem  which  is  the  adjoint  of  the  original  boundary  value  problem. 

In  terms  of  the  original  and  adjoint  problems,  we  will  then  be  able 
to  state  a  functional,  the  stationary  value  of  which,  in  the  sense  of 
the  calculus  of  variations,  yields  the  differential  equations  and 
boundary  conditions  for  both  problems.  In  [19]  it  was  stated  that 
there  is  some  evidence  that  this  so-called  adjoint  variational  prin¬ 
ciple  leads  to  slightly  better  results  than  does  the  Galerkin  method 
and  that  the  speed  of  convergence  might  also  be  better.  As  is  perhaps 
to  be  expected,  the  somewhat  more  accurate  results  obtained  with  the 
adjoint  variational  principle  are  acquired  at  the  expense  of  increased 
complexity  and  labor. 

The  importance  of  the  adjoint  problem  in  nonself-sdjoint  eigen¬ 
value  problems  was  discussed  by  Roberts  [21],  Ballio  (??]  -  [23] 
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applied  the  concept  of  the  adjoint  problem  to  a  couple  of  nonconserva¬ 
tive  stability  problems,  including  Beck's  problem.  Nemat-Nassar  and 
Herrmann  [24]  also  determined  the  boundary  value  problem  adjoint  to 
beck's  problem,  and  Prasad  and  Herrmann  [25]  described  the  usefulness 
of  tlie  adjoint  problem  in  solving  nonconservative  stability  problems. 
The  method  described  in  [25]  has  been  generalized  by  Anderson  and 
Walter  [26]  and  .Anderson  [27]  to  include  the  effects  of  internal  and 
external  damping  in  nonconservative  stability  problems. 

It  was  shown  in  [24]  -  [26]  that  the  eigenvalue  problem  that  is 
adjoint  to  (6)  and  (7)  is 

v^(x)  +  Q  v**(x)  -  a)^v(x)  a  0,  o<x<l, 

(29) 

v(o)  =  V  (o)  =  v'  (1)  +  Q  v(l)  =  v'**(l)  +  Q  v*(l)  =  0, 

(30) 

which  is  often  called  Peut's  problem,  .Moreover,  it  was  shown  in 
[25],  [26],  and  [28]  that  both  the  original  problem  (6),  (7)  and 
the  adjoint  problem  (29),  (30)  can  be  generated  from  the  variational 
principle 

6(  /^[y"v*'  -  Qy  v'  -  u^yv]dx  Qy  (l)v(l)}  =  0, 

0  (31) 

where  now  we  are  dealing  with  the  problim  of  finding  a  stationary 

value  tor  the  given  functional.  We  shall  now  use  (31)  in  place  of 

(9)  for  the  purpose  of  determining  approximate  values  for  the  critical 

load  in  Beck's  problem.  It  is  essential  to  keep  in  mind  that  the 

eigenvalues  of  the  adjoint  problem  are  identical  to  the  eigenvalues 

of  the  original  problem. 
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To  obtain  a  suitable  set  of  coordinate  functions  for  the  adjoint 
problem,  let  us  assume  that 
^  N 

V(x)  =  v(x)  =  I  b  V  (X). 

n=l  "  (321 

Moreover,  since  the  geomecric  boundary  condltioiis  in  (30)  are  the  saine 

as  those  in  (7),  we  shall  assume  that  the  Vj^(x)’s  are  of  the  form  shown 

in  (11),  i.e., 


Vn(x)  =  ^  n  >  1. 


If  we  now  require  (33)  to  satisfy  the  two  boundary  conditions  at  x  =  1 
that  are  stated  in  (30),  we  can  show,  following  the  method  applied 
above  to  (11)  and  (12),  that  the  coefficients  Aj^,  B^,  can  be  selected 
such  that  (33)  assumes  the  form 


v„(x)  =.  v"*‘  -  *  V"*'- 


where 


Oy,  =  ♦  2n(n+2)  Q  +  (n+l)(n*2)^(n+3). 


2n(n+l)  Q  +  n(n+l)  (n'«-2)  (n+3) , 


Yjj  =  +  2(n^-l)  Q  +  n(n+l)^(n+2). 


If  we  now  si'bstitute  (10)  and  (32)  into  (31),  we  obtain 
N  N 

'  I  I  W  “nk  "  0. 

n^l  k=l 
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whei*e 


WVj^' (x)-Q  y^(x)vj^(x)-w^y^(x)Vj.  (x)]dx  +  Q  y^CDv^Cl). 
o 

(36) 

Performing  the  variation  in  (35),  we  find 

\  ‘„k  '"k  *  ‘kn  ‘’k)  =■ 

n=l  k=l 

from  which,  in  the  usual  way,  it  follows  that 


1  ‘nk  '  “•  1  '’n  ’  »• 


n=l 


n=l 


(37) 


The  frequency  equations  are  then 

det  (I^k^  "  ^  (^kn^  " 


(38) 


But  these  equations  are  equivalent  because  the  value  of  a  determinant 
remains  unchanged  upon  interchanging  its  rows  and  columns.  This  is 
simply  a  manifestation  of  the  fact  that  the  eigenvalues  of  the  original 
and  adjoint  problems  are  identical. 

Using  tne  coordinate  functions  in  (13)  and  (34),  which  satisfy  all 
the  geometric  and  natural  boundary  conditions  in  (7)  and  (30),  we  may 
use  either  determinant  in  (38)  to  evaluate  the  critical  load.  The 


value  reported  in  [25]  and  [28]  is  for  N  =  2 


Qcr  =  1.97nr% 


(39) 


somewhere  lower  than  the  exact  value  given  in  (8). 

Let  us  go  one  step  further  without  consideration  of  (36)  and  (38), 
Recalling  that  in  the  application  of  the  Ritz  method  with  a  functional 
of  the  type  given  in  (31),  we  are  not,  in  general,  required  to  select 
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coordinate  functions  which  satisfy  both  the  geometric  and  natural 
boundary  conditions  of  the  given  eigenvalue  problem,  indeed,  we  are 
really  compelled  to  satisfy  only  the  geometric  boundar)'  conditions, 
although  this  usually  means  that  the  speed  of  convergence  of  the 
method  is  slowed  somewhat  depending  upon  the  nature  of  the  problem 
under  consideration.  Let  us  now  assume  that  the  functions  YjjC^)  and 
Vj.(x)  in  C36)  satisfy  only  the  geometric  boundary  conditions 


YnCo)  =  yn(o)  =  v^Co)  =  v^(o)  =  0. 


(40) 


In  particular,  we  are  assuming  now  that  y^^Cx)  and  v^^Cx)  ?re  not 
necessarily  the  functions  defined  above  in  (13)  and  (34). 

Integrating  by  parts  and  using  (40),  we  can  easily  demonstrate 

that 


/V^’(x)Vj^’(x)dx  =  /V^  (x)Vj^(x)dx+yj^'COVj[(l)-y^"(l)v^(l), 


=  -  /Vn'(x)Vj^(x)d.4  *  /n(^)\(l)» 


o’'"  0 

so  that  (36)  becomes 


Therefore,  the  first  relationship  in  (34)  can  be  written  as 
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or,  equivalently 


Ay®^Cx)  +  Q  y"(x)  -  i/y(x)]v,ix)dx  + 


+  y"  (l)vj^(l)-y' "  (l)v^(l)  =  0. 


(41) 


liq.  (41)  bears  a  very  strong  resemblance  to  the  so-called  semi-extrnded 
Galerkin  equations  for  nonself-ad joint  eigenvalue  problems  as  discusi?ai 
by  Leiphol;  [2],  Note  that  here,  however,  the  weight  functions  are  tho 
coordinate  functions  of  the  adjoint  problem  V|^(x)  rather  than  the 
coordinate  functions  of  the  original  eigenvalue  problem.  The  presence 
of  the  terms  y'*(l)vj!(l)  -  y*  *  *  (l)Vj^(l)  in  (41)  serves  to  compensate 
for  the  fact  that  tfie  coordinate  functions  yjj(x)  do  not  satisfy  the 
natural  boundary  conditions  y* *  (1 )  =  y***{l)  =  0  of  the  original 
eigenvalue  problem,  lienee,  (41)  may  be  considered  a  generalization 
of  the  results  presented  in  [2], 

If  the  coordinate  functions  y„(x)  satisfy  the  natural  boundary 
conditions  of  the  original  problem,  then  (41)  reduces  to 


[y^(x)  Qy"(x)  -  u‘^y(x)]v  (x)dx  =  0, 


2'^'. 


(d2) 


which  is  analogous  to  the  classical  fonri  of  the  Galerkin  equations. 

To  this  point  we  have  assumed  that  t’-.e  weight  functions  V|^(x)  satisfy 
only  the  geometric  boundary  conditions  of  tne  adjoint  prob<em.  The 
question  of  whether  the  procedure  embodied  in  (42)  converges  requires 
further  study.  If  this  procedure  does  converge,  it  is  very  probable 
that  it  will  converge  even  more  rapidly  if  the  Vj^(x)'s  are  also  *e- 
quirod  to  satisfy  the  natural  boundary  conditions  of  the  adjoint  problem. 
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If  and  v^(x)  are  defined  by  (13)  and  (34),  then  (42)  is 

equivalent  to  (37).  In  general,  the  numerical  work  is  less  complicated 
in  (42)  then  in  (37)  with  I^^j^  given  by  (36), 

THE  COMPLEMENTARY  ENERCY  AND  ADJOINT  VARIATIONAL  METHODS 


Now  let 

M(x)  =  -  y"  (X),  N(x)  =  -  v"  (x) 

so  that  (31)  can  be  expressed  as 


(43) 


6{/^ [M(x)N(x)  -  0.  y*(x)v*(x)  -  w“y(x)v(x)]dx  + 
0 

Q  y’(l)v(l)}  =  0. 


(44) 


The  "bending  moments"  M(x)  and  N(x)  are  subject  to  the  auxiliary 
conditions 


m"(x)  =  Q  y"(x)  -  wV(x), 

M(l)  =m'(1)  =  0, 

n”  (x)  =  Q  v"  (x)  -  u^v(x), 

N(l)  =Qv(1),  n’(1)  =  Q  v'(1). 


(45) 


(46) 


Because  we  have  introduced  the  adjoint  problem  and  consequently  a 
variational  principle,  we  may  view  (44)  as  a  replacement  for  f21)  with 
(45)  and  (46)  being  introduced  in  analogy  with  (??)-(23).  References 
[10l-[i2)  have  again  supplied  the  motivation  for  this  type  of  analysis. 
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We  shall  again  assume  that  the  coordinati  functions 
Vj^(x)  are  those  given  in  (13)  and  (34).  Since  (24)-(27)  are  still 
applicable,  we  need  only  derive  an  expression  for  N^(x)  before  pro¬ 
ceeding  to  the  numerical  work.  To  this  end,  we  find,  after  integrat¬ 
ing  the  first  equation  in  (46)  and  imposing  the  boundary  conditions 
in  (46), 

N(x)  Q  v(x)  -  (x-t)v(t)dt. 

1 


If  we  insert  (32)  into  this  last  result,  we  obtain 
N 

N(x)  =  I  bn  Mn(x). 

n=l  (47) 

whore 

Nn(x)  =  Q  Vn(x)  -  f  (x-t)Vn(t)dt 

1 


=  Q  v„(x)  -  w\(x), 

(48) 

with 


2(q^+3(n2*n-4)Q+2(n-^l)  (n+2)  (n+3)  (3n+10)] 

t'n(x)  =  /  f:--t)Vj^(t)dt  =  — - - 

1  (n+3) (n+4) (n+5) 


2x[Q^-t-5(n-l)  (n-^2)Q-^6(n-H)  (n*2)"  (n»3)l  + 

(n+2)  (n+3)  (n+4) 


(  (n+2)  (n+3)  (n+3)(n+4) 


(n+4)  (n+5) 
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Substitution  of  (26),  (47),  (10),  and  (32)  into  (44)  yields 


N  N 

^  I  I  ®n  '^nk  "" 

n=l  k=l 


(49) 


where 

Jnk  -  /^[Mn(x)‘\(X)-Q  yj^(x)v^(x)-wVj^(x)Vj^(x)]dx+Q  yj^(l)Vj^(l) 
0 


=  /^M^(x)N^(x)  ♦  [Q  y^*(x)-u^.^(x)]Vj^(x)}dx. 


In  view  of  (13),  (27),  (34),  and  (48),  the  coefficients  may  now  be 
considered  completely  known.  Proceeding  with  (49)  as  we  did  with  (35), 
we  find 

N 

y  J  ,  =  0, 

n  nk  * 
n=l 


and  consequently 

det  (J  ,  )  =  0, 

nk 

Numerical  calculations  with  (50)  yield,  for  N  =  2, 


(50) 


=  2.061^^ 


(51) 


the  most  accurate  of  all  the  approximate  calculations  performed  here. 

To  summarize  the  numerical  results  obtained  by  the  various  methods 
described  above,  the  values  of  are  now  tabulated  in  Table  I  along 
with  the  percentage  of  error  relative  to  Beck's  exact  value. 
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TABLE  I 


COMPARISON  OF  THE  VARIOUS  METHODS 


EQUATION  NUMBER 

(8)  Exact 

.. 

(16) 

7,98 

(28) 

3.15 

(39) 

1.971 

-2.96 

(51) 

2.061 

1.48 

Examination  of  Table  I  reveals  that  the  Ritz  procedure  based  upon 
the  variational  equation  in  (9)  yields  the  poorest  result,  which  is 
almost  8o  {greater  than  the  exact  value  for  Using  the  variant  (21) 

of  (9),  i.e.,  the  method  of  complementary  energy,  we  obtain  a  much  bet¬ 
ter  approximation,  namely  (28),  which  is  about  3%  greater  than  the  exact 
value.  A  slightly  better  approximation,  (39),  was  obtained  by  means  of 
the  Ritz  method  based  upon  the  adjoint  variational  principle  in  (31), 
However,  the  best  approximation  of  all,  (SI),  was  obtained  through  a 
scheme  which  combined  the  methods  of  the  adjoint  variational  principle 
and  the  complementary  energy.  In  this  later  case,  the  approximate  value 
of  Qj.j.  was  only  about  1  l/2c  greater  than  the  exact  value.  It  must  be 
admitted  that  the  improved  value  given  in  (51)  was  obtained  at  the 
expense  of  (i)  introducing  the  adjoint  problem  (29)  -  (30)  and  the 
adjoint  variational  principle  (31),  (ii)  selecting  a  set  of  coordinate 
functions  for  the  adjoint  problem  (34),  and  (iii)  deriving  expressions 
for  the  "bending  moments"  integrating  (45)  and  (46), 

It  might  be  argued  from  the  standpoint  of  the  labor  involved  that  it  is 
more  efficient  to  use  the  ordinary  Galerkin  procedure  (17)  with  N>3. 


flAUGER'S  PROBLEMS 


In  the  previous  paragraphs  we  discussed  Beck's  problem  for  which 
it  is  possible  to  obtain  the  exact  value  of  the  critical  load.  Using 
certain  methods  of  approximation,  we  calculated  approximate  values 
for  and  compared  these  results  with  the  known  exact  value.  In 
many  nonconservative  stability  problems  it  is  a  very  difficult  task 
to  solve  the  differential  equations  of  motion  exactly,  and  consequently 
the  value  of  the  critical  load  parameter  can  be  found  only  with  consid¬ 
erable  effort.  This  is  particularly  true  for  differential  equations 
that  have  variable  coefficients.  In  such  cases  the  methods  described 
in  the  preceding  paragraphs  offer  rather  efficacious  tools  for  solving 
nonconservative  stability  problems.  Hauger  [29]  considered  problems 
of  this  type,  namely,  he  computed  approximately  the  critical  loads  for 
beams  subjected  to  (i)  uniformly  distributed  and  (ii)  linearly  distrib¬ 
uted  tangential  compressive  loads. 

We  shall  concern  ourselves  here  with  two  problems  studied  by  Hauger 
in  which  a  linearly  distributed  tangential  load  is  acting  along  the 
length  of  the  beam.  If  the  boundary  conditions  for  the  beam  are  either 
one  end  clamped  -  one  end  free  or  one  end  clamped  -  one  end  hinged,  then 
the  kinetic  method  must  be  used  to  determine  the  critical  load  because 
these  systems  become  unstable  through  flutter.  The  equation  of  motion 
is 


3^w  1  ..  -.2 

3x^  2  9t^ 


(52) 
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Merc  denotes  the  intensity  of  the  tangential  load  and  the  remaining 
quantities  are  identical  with  those  defined  in  (1).  If  we  introduce 
the  same  length  and  time  parameters  that  we  used  to  derive  (3),  then 
(52)  can  be  expressed  as 
,4 


where 


2  a^w 

— ^  ♦  Q(l-x)'^  —  —  =  0, 

ax^  3x  at 


^0^ 

2E1 


0<X<1,  t>0, 


(53) 


Setting  w(x,t)  =  y(x)  e^^^,  we  reduce  (53)  to  the  following  ordinary 
differential  equation  with  variable  coefficients: 


y^(x)  +  Q(l~x)^  y’*(x)  -  uV(x)  =  0,  o<x<l. 


For  a  clamped-free  beam  the  boundary  conditions  are 


y(o)  =  y'(o)  =  y''(l)  =  y”'(l)  =  0 


and  for  a  clamped-hinged  beam  they  are 


(54) 


(55) 


y(o)  =  y' (o)  =  y(l)  =  y’'(l)  =  0. 


(56) 


Henceforth,  we  shall  refer  to  (54)  subject  to  (55)  as  Hauger's  first 
problem  or  Problem  I  and  to  (54)  subject  to  (56)  as  Hauger's  second 
problem  or  Problem  II, 

because  of  the  difficulty  of  solving  these  two  boundary  value 
problems  exactly,  Haugo^  resorted  to  Galerkin's  procedure  in  its 
classical  form,  employing  clamped-free  and  clamped-hinged  beam  modes 
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(obtained  from  (54)  -  (56)  with  Q  =  0)  for  the  coordinate  functions. 
According  to  Leipholz  investigations  [4]  -  [7],  this  procedure  should 
converge  to  the  exact  value  of  the  critical  load  intensity  pro¬ 
vided  that  a  sufficient  number  of  terms  is  taken  in  the  approximation 
for  y(x).  liauger  used  only  two-teimi  approximations  and  calculated  the 
following  values  for  the  critical  load  intensity  parameter 


2  j 

Qcr  =  8.0171  and  Q  =  20.3877^, 


for  Problems  I  and  II,  respectively. 

Il  is  our  intention  here  to  reconsider  these  problems  and  to 
demonstrate  by  numerical  calculation  that  these  values  are,  respec¬ 
tively,  4.98^.  and  28.3%  greater  than  the  correct  values.  In  view 
of  the  outcome  of  the  calculations  tabulated  in  Table  I  for  Beck's 
nroblem,  we  shall  elect  to  introduce  the  appropriate  form  of  the 
adjoint  variational  principle  for  Hauger's  problem.  The  numerical 
calculations  to  be  presented  here  will  be  based  on  the  methods  of 
Halerkin  and  complementary  energy  in  conjunction  with  the  adjoint 
variational  principle. 


THE  GALERKIN  PROCEDURE  AND  THE  ADJOINT  VARIATIONAL  PRINCIPLE 


It  was  shown  in  [27]  that  the  differential  equation  that  is  the 
adjoint  of  (54)  is 

v^(xj  Q  [  (l-x)^v(x)] '  *  -  a)^v(x)  =  0,  o<x<l, 

(58] 

while  the  boundary  conditions  for  the  adjoint  problem  are  identical 
to  (55)  and  (56)  for  Problems  I  and  II. 
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Furthermore,  it  was  shown  that  these  two  boundary  value  problems  may 
be  generated  from  the  variational  principle 


5/‘{y' ’  (x)v* '  (x)  -  Q  y*  (x)  ( (l-x)^vCx)] '-a)^y(x)v(x)  }dx  =  0. 


We  next  assume  that 


y(x)  y(x)  =  I  a  y  (x),  w(x)  =  v(x)  =  I  b  v  (x), 
n=l  n*l  " 


where,  for  the  present,  the  coordinate  functions  Xj^Cx),  Vj^(x)  are 
assumed  to  satisfy  only  the  geometric  boundary  conditions  of  a  given 
problem,  lluis,  for  Hauger's  first  problem,  we  assume  that 


YnCo)  =  XnCo)  = 


whereas  for  Hauger's  second  problem 


/nCo)  =  y^(o)  =  y^(l)  =  v^{o)  =  Vjj(o)  =  v^^Cl)  =  0. 


Substitution  of  (60)  into  (59)  yields 


N  N 


&l  I  a„b^  /^{y’’(x)v’'(x)-Q  y’(x)[(l-x)2v  (X))’ 
n=l  k=l  0 


-  “  y„ (x)v.  (x) }dx  =  0. 

K 


Integrating  by  parts,  we  can  easily  show  that 


/  y^^  (x)Vj^(x)dx  +  y^  -  y^’  (l)v^(l). 


(Problem  I) 


/V^’  fx)Vj^'  (x)dx  =1 


/V^(x)v  (x)dx  +  y’’(l)v'(l), 

.0  "  HR  (Problem  II) 


y' (x) f  (l-x)^v  (x)]  dx  = 
o  " 


-  (l-x)V'(x)v  (x)dx, 

o  n  K 


where  we  have  imposed  the  geometric  boundary  conditions  (61)  and  (62), 
Inserting  these  last  two  results  into  (63),  we  obtain 


N  N 


i  ‘<[y_(x)]Vj.(x)dx  ♦  y’'(l)v’(l)-y;'*(l)v^(l)}  =  0 

n=l  k=l  o 

(64) 


M  N  1  , 

^  I  I  a^b^  {/  L[yj^(x)]v^(x)dx  +  y^'(l)v^(l)}  =  0 
n=l  k=sl  0  (65) 


for  Problems  I  and  II,  respectively,  where 


bfyj^(x)]  =  *  Q(i-x)  y^^’ (x)  -  w^yj^(x). 


Proceeding  with  (64)  and  (65)  in  the  usual  fashion,  we  arrive  at 


I  (/^  b[y^(x)]v^(x)dx  +  y^’(Uv^(l)  -  y^' (l)v^(l)}  =  0 


n-trn  -  . ’k 


'n 


2'.) 


I  {/  L[y^(x)]v^(x)dx  ♦  =  0 

n=l  o 


or  equivalently 
rl  . 


/  L[yCx)lv^(x)dx  ♦  y’  Cl)vj^(l)  -  y  (l)Vj^(l)  «  0, 


Uy(x)]v^(x)dx  >  y”a)v,^(l)  =  0 


for  Problems  I  and  II,  respectively,  Eqs.  (66)  and  (67)  are  in  the 
"semi-extended"  form  of  (5alerkin‘s  procedure  as  described  in  |2] 
except  that  the  weight  functions,  Vj^(x),  here  arc  the  coordinate 
functions  of  the  adioint  problem  rather  than  of  the  original  problem 
as  is  normally  the  case. 

Because  the  speed  of  convergence  is  faster,  in  general,  when  the 
coordinate  functions  are  selected  so  as  to  satisfy  all  the  geometric 
and  natural  boundary  conditions  of  a  given  problem,  we  shall  elect  to 
meet  this  requirement,  thereby  reducing  (66)  and  (67)  for  both 
Problems  I  and  II  to 

L[y(x)]v.  (x)dx  =  0. 


Moreover,  since  the  boundary  conditions  of  the  original  and  adjoint 
problems  are  identical,  we  shall  select  the  y^^’s  and  Vjj’s  to  be 
identical.  Then  (68)  becomes 

L[y(x)]yk(x)dx  =  0,  (69: 

0 

which  is  the  Galerkin  procedure  in  its  familiar  form.  Our  numerical 
calculations  will  be  based  upon  (69). 
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ror  Problem  I,  the  clamped- free  beam,  the  coordinate  functions  in 


(13)  are  again  suitable.  Numerical  caJculationr  were  performed  with 
(69)  and  (10)  for  N  =  2, 3,..., 7,  and  the  results  are  shown  .in  Tvable  II. 
TABLE  II  -  THE  CRITICAL  LOAD  INTENSITY 
FOR  iiAUGER’S  FIRST  PROBLEM 


N 

2 

7.5776 

3 

7.7830 

4 

7.6446 

5 

7.6310 

6 

7.6318 

7 

7.6314 

For  Problem  II,  the  clamped-hinged  beam,  the  set  of  polynomial 


coordinate  functions 

y  (x)  =  (n+2)x"*^-  (2n+3)x"'*'^+  (n+])x"'’^,  n>l 
"  “  (70) 


satisfies  all  the  geometric  and  natural  boundary  conditions  in  (56). 


Again  using  (69)  and  (10),  we  obtain  the  results  exhibited  in  Table  III 


TABLE  III  -  THE  CRITICAL  LOAD  INTENSITY 


FOR  HAUGER’S  SECOND  PROBLEM 


THE  COMPLEMENTARY  ENERGY  METHOD 


Before  commenting  on  the  numerical  results  presented  in  Tables  II 
and  III  in  relation  to  Hauger's  values  in  (57),  we  shall  re-solve  these 
problems  using  now  the  method  of  complementary  energy.  Introducing 
again  (43) ,  we  may  express  (59)  as 


6/\m(x)N(x)-Q  y*  (x)  [(l-x)^v(x)]*-u)^y(x)v(x)}dx  =  0, 

0  (71) 

with  the  auxiliary  conditions 

m"(x)  =  Q(l-x)2y"  (x)  -  (ij^y(x), 

(72) 

n"(x)  =  Q((l-x)^v(x)] "  -  u)2v(x), 

subject  to 

M(1)  =  m’(1)  =  N(l)  =  N'(i)  =  0 

(73) 


for  Problem  I  and 


M(l)  =  N(l)  =  0 


(74) 


for  Problem  II,  which  are  consistent  with  (54)  -  (56),  (58). 

For  Hauger's  first  problem  we  shall  assume  that 

N  N 

y(x)  =  I  v(x)  =  I  b  yj^(x), 

n=l  n=l  (75) 

where  the  yj^(x)'s,  which  satisfy  all  the  boundary  conditions  in  (55), 

are  given  in  (13),  By  virtue  of  (72)  and  (73),  we  may  also  write 


N  N 

M(x)  =  I  a^M  (x),  N(x)  =  I  bj^Njj(x), 
n=l  n=l 

(76) 
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where 


=  Q(l-x)2y^'(x)-toVn(x),  Mj^(l)  =  M^^Cl)  =  o. 


(77) 


M  =  Qi(l-x)  yn(x)]'’-u:2y^(x),  N^(l)  =  Nj^(l)  =  0. 
Integration  of  (77)  yields 

Mn(x)  =  Q[(l-x)2y^(x)  +  2(x+2)r„(x)  -  6Sj^(x)]  -  w^gj^(x), 
.\(x)  =  ()(l-x)2y^(x)  -  i/g^(x), 
where 


(78) 

(7P) 


1 


_6_  ^ 

n+4  ^  2 

^  n+3  -  2nx  + 

n(n+l) 

n+4 

2  (3n-fl0)  f 

(n-*'4)  (n+5)  2  \ 

2n  (n+3) 

n+4 

2f3ni-101 

6x  1 

n+4 

(n+4) (n+5) 

n+4  2 

(76)  into 

(71),  we  find 

N  N 

*  I  X\\  'nk  '  0. 


n=l  k=l 


where 


^nk  ■  ^  =  Q(l-x)^y„' (x)yj_  (x)  -  w^y^(x)yj^(x)]dx. 


This  leads  to 
N 

I  \  ^nk  =  0- 


n=l 


(80) 


Tne  numerical  results  obtained  from  (80)  are  shown  in  Table  IV. 
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TABLE  IV  -  THE  CRITICAL  LOAD  INTENSITY 
FOR  HAUGER’S  FIRST  PROBLEM 


N 

Q 

xr 

2 

7.573 

3 

7.653 

4 

7.631 

5 

_ 

7.631 

For  Hauger's  second  problem  we  may  again  adopt  (75),  however  now 
the  coordinate  functions  will  be  those  which  appear  in  (70).  Then 
with  (70)  apd  (75)  integration  of  (72)  subject  to  (74)  leads  to 


N 

M(x)  *  aQ(x.-i)  +  I  anMj^(x), 
n=l 

M 

N(x)  =  bo(x-i;i  +  I  b„Nn(x), 
n=l 


(81) 


where 

a^  -  m'(1),  bg  =  N'd). 

and  M^(x)  and  Nj^(x)  are  given  by  (78)  and  (79),  respectively,  with 


r„(x)  =  /  y„Ct)dt  = 


-3 


(n+3)  (n 


_.x"*2(i.£2ia.x.JS!iax2) 

+4)  \  n+3  n+4  ' 


Sn(x)  =  /  t  y  (t)dt 
1 


■(5n+7) 


(n+3)  (n+4)  (n+5) 


+  X 


n+3^n+2 
ln+3 


2n+3  t;+ 1 

- r  X  ♦  — r  X  j, 

n+4  n+S  / 
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Sn'*’? 


3x 


(n+5)  (n+4)(n+5)  (n+3)(nM) 


+  {  -i-  -  +  (n^l)x^ 

(n+3)(n+4)  (,->+4)(n*5) 


Inserting  (75)  aii<l  (81)  into  the  variational  principle  (71),  we 


find 


1  •^'  N  N  N 

-  l>0  I  Vn  *  I  I  V’klnk'-O. 

n=l  n=l  k=l 

(82) 


n=l 


whore 


‘nk  ""  /  ■*■  <Kl~x)^yj^'(x)yj.(x)-a)“y^(x)y  (x))dx, 


(l-x)Mj^(x)dx,  (l-x)N„(x)dx. 


hut  (82)  leads  to 


N  S 


T-o-  IVn  “'’o*  !  K>„k-Vk  ^l-k 

*  k—  I  n—  1 


N  N 


"  I  i^k^nk-Vn  ^^n  =  0. 

n=l  n=l  k-1 

(83) 

In  the  usual  fashion,  we  obtain  from  (83) 


‘^0  '  ^  Vn  =  0.  I  ^  Ink  “  “o^k  = 

n=l  n=l 


(84) 


plus  two  additional  equations  which  can  easily  be  shown  to  be  equivalent 


to  (84).  Eliminating  a^  between  the  two  equations  in  (84),  we  arrive  at 
N 

n=l  "  (85) 

The  numerical  calculations  obtained  from  (85)  are  shown  in  Table  V. 

TABLE  V  -  THE  CRITICAL  LOAD  INTENSITY 
FOR  HAUGER'S  SECOND  PROBLEM 


N 

2 

16.522 

3 

16.012 

4 

15.929 

5 

15.884 

Therefore,  in  view  of  Tables  II  -  V,  we  conclude  that  accurate 
values  for  the  critical  load  intensity  parameter,  are 

2  2 
Q  =  7.631TI  and  =  15.88ti^ 

cr  cr 

for  Ilauger's  first  and  second  problems,  respectively.  These  results 
are  4.98’o  and  28. 3o  lower  than  the  corresponding  numerical  values 
reported  by  Hauger  (see  Eq,  57)  who  used  a  two-term  Galerkin  approxi¬ 
mation  with  clamped-free  and  clamped-hinged  beam  modes  as  coordinate 
functions.  Here  for  coordinate  functions  we  have  used  polynomials 
that  satisfy  all  the  geometric  and  natural  boundary  conditions  of  the 
jspective  problems  in  (i)  the  Galerkin  procedure  and  (ii)  the  method 
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of  complementary  energy  in  conjunction  with  the  so-called  adjoint 
variational  principle  to  calculate  approximate  values  for  The 

anpli cation  of  these  processes  was  continued  by  successively  increas¬ 
ing  the  number  of  terms  included  in  the  approximation  functions  until 
the  numerical  values  appeared  to  converge  to  the  correct  value  for 
Q^j..  In  both  cases  it  appears  that  the  method  of  complementary  energy 
converged  at  a  somewhat  faster  rate  than  did  the  Halerkin  procedure. 

Hie  r.cthod  of  complementary  energy  requires  extra  effort  to  apply 
relative  to  the  Galcrkin  method  since  the  "bending  moments"  must  be 
obtained  by  integrating  an  elementary  second  order  linear  differential 
equation.  The  Galcrkin  procedure  offers  the  advantage  that  the  numeri¬ 
cal  calculations  may  be  set  up  in  a  rather  straightforward  fashion. 

\t  least  in  the  numerical  computations  reported  hero,  the  application 
of  the  Galcrkin  procedure  offered  no  computational  hardships  as  the 
value  of  the  integer  N  increased  because  all  the  integrals  involved 
were  evaluated  very  efficiently  on  the  computer  by  means  of  a  double 
precision  subroutine  based  upon  the  Gaussian  quadrature  technique 
which  evaluates  integrals  of  polynomials  exactly.  .Vnother  sul)routine 
was  used  to  evaluate  the  various  determinants  that  arose,  and  a  root 
finding  scheme  determined  the  value  of  the  load  intensity  parameter 
for  a  given  value  of  the  frequency  of  vibration  of  the  system. 
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